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Abstract. We present an analysis of the magnetic- field fluctuations in the 
magnetoionic medium in front of the radio galaxy 3C 31 derived from rotation- 
measure (RM) fits to high-resolution polarization images. We first show that the 
Faraday rotation must be due primarily to a foreground medium. We determine 
the RM structure functions for different parts of the source and infer that the 
simplest form for the power spectrum is a power law with a high-frequency 
cutoff. We also present three-dimensional simulations of RM produced by a 
tangled magnetic field in the hot plasma surrounding 3C 31, and show that the 
observed RM distribution is consistent with a spherical plasma distribution in 
which the radio source has produced a cavity. 



1. Observed Faraday rotation and depolarization 

Our knowledge of the structure and origin of magnetic fields in elliptical galax- 
ies, groups and clusters is still rudimentary, but Faraday rotation of linearly- 
polarized radio emission can be used to probe the fields in ionized foreground 
gas. Here, we present an analysis of the magnetic-field fluctuations in the mag- 
netoionic medium in front of the FRI radio galaxy 3C 31 (z = 0.0169) derived 
from rotation-measure (RM) fits to high-resolution polarization images. 

Our analysis is based on VLA observations at 6 frequencies in the range 
1.4 - 8.4 GHz with resolutions of 5.5 and 1.5arcsec FWHM (Figs Q] and [2]) . We 
show images of normalized polarization gradient p'(0)/p(0) from fits to p(X 2 ) = 
p(0) + p'(0)A 2 , where p{\ 2 ) is the degree of polarization at wavelength A and 
a prime denotes differentiation with respect to A 2 , together with RM images 
derived from fits to x(X 2 ) = x(0) + RMA 2 at 4 - 6 wavelengths, where x is the 
E-vector position angle. 

The residual depolarization at 1.5arcsec resolution is very small and the 
rotation is accurately proportional to A 2 , indicating almost completely resolved 
foreground rotation. There is a large asymmetry across the nucleus: the lobe 
with the brighter jet shows a much smaller RM fluctuation amplitude than 
the cou nter-jet lobe on all scales, qualitatively as expected from relativistic jet 
models (|Laind[i~988l ). 
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Figure 1. Images at 5.5arcsec resolution, (a) Mean of four total- intensity 
images at frequencies from 1.4 - 1.7 GHz. (b) Normalized polarization gradi- 
ent j>'(0)/p(0), in m -2 . (c) RM from a fit to 4 images at frequencies between 
1.4 and 1.7 GHz, in radm -2 . (d) As (c), but for 5 frequencies from 1.4 - 5.0 
GHz. 
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Figure 2. Images at 1.5 arcsec resolution, (a) Mean of 4 total-intensity im- 
ages at frequencies from 1.4 - 1.7 GHz. (b) Normalized polarization gradient 
p'(Q)/p(Q), in m~ 2 . (c) RM in radm -2 from a fit to 6 frequencies between 1.4 
and 8.4 GHz using the PACERMAN algorithm (|Dolag. Vogt fc En61inll2005f) . 
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Figure 3. (a) - (f): plots of the RM structure function for the regions 
indicated in Figs Q] and [2] The curves are the predictions for the cut-off 
power law power spectrum described in the text, including the effects of the 
convolving beam. 



2. Structure functions and power spectra 



We quantify the spatial variation of RM using the structure function S(r) = 
([RM(r + r') — RM(r')] 2 ), where r and r' are separation vectors in the plane 
of the sky and () denotes an average over r'. Our data are consistent with 
isotropy, so we plot 5 as a function of scalar separation r for five regions of 
3C 31 in Fig. [3l In order to interpret these results, we start from a model 
RM power spectrum C(k), where k is the wave-number. Rs Hankel transform 
is the RM autocorrelation function C(r), and the structure function is S(r) = 
2[C (r) — C (0)] . A novel feature of our analysis is that we include the effects of the 
observing beam explicitly (this can be done straightforwardly provided that the 
rotation across the beam is small, as is the case for our observations). Our data 
are consistent with a power spectrum which has the same form everywhere, but 
varying normalization. We have found two acceptable functional forms for the 
power spectrum: C(k) oc /c -2 ' 35 with a high-frequency cut-off at k = 0.5 arcsec -1 
(a scale of 12 arcsec or 4kpc) and a broken power-law form with C(k) oc /c~ n / 3 
(as expected for Kolmogorov turbulence) for k > 0.13 arcsec -1 and oc k at 
larger scales. The predicted structure functions are very similar if the effects 
of the beam are included. This may explain why earlier studies have come to 
different conclusion s regarding the index of the power spectrum on the basis of 
RM analysis alone dVogt fc EnBlinll2003l . 120051 : iMurgia et aDl2004l : iGovoni et af 



2006). The cut-off power- law model with C{k) oc k 2 ' 35 predicts significantly less 



residual depolarization, however, in better agreement with our data for 3C31. 



4 




-50 



I 1 1 1 I 



— : — ~ — r -1 — i — 1 — r~ 

(b) Mean HM 
Observed 



J I I I I I I I I ■ ■ 1 ! I I I I L 



id 



3(« 



1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


- 1 — 1 1— 1 


Ot served - 


1-4 1 — 1 








~ •— 1 

"I , . . 1 . , , 


i | ^ ^ ' — 1 " 
, . . i . . . i ." 




_1 I I I 1_ 



so 



,,,,,,, 

■ 


1 

(e) Mean RM _ 
No cavity . 

■ 







1 — 1 
1 — 1 


_^fc=T± 

■ 




... i ... i . 



iv:: 



fit! 



in 



1 

1 — 1 


, , , 1 , 

No cavity - 


S \ — i 

■l , . . 1 , , , 


, . . i . . . 1 .- 




I 1 1 1 I 



(h} Mean RM 
Cavity 



-r=t?*t 



l ■ . ■■ l 



J I I 1 I I I L 



J 1 1 1 1 1 1 1 




'_ 1 — 1 


(!) "b, 


M 


Cavity 


- I_l i— i 1— 


— i : 


- 

'' 1 ' ' 1 ' 1 1 


i— i— ■ 

, , , i . . . 1 ,' 



-400 -$ty 200 400 
Distance Irom nucleus / arcsec 



-400 -aoo o aoo 40c 

Distance from nucleus / arcsec 



Figure 4. Comparison of RM data at a resolution of 5.5 arcsec FWHM 
with three-dimensional simulations. Panels (a) - (c): data; (d) - (f) simula- 
tion with spherically symmetrical gas distribution (central rms magnetic field 
strength (S 2 ) 1 / 2 = 0.21 nT); (g) - (i): simulation with cavity, (a), (d) and 
(g) show RM images, (b), (e) and (h) show mean values over boxes of length 
65 arcsec along the jet axis, (f) and (i) the corresponding rms variations. The 
curve in panel (f) is the variation predicted by a single-scale model. 
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3. Three-dimensional simulations 

The easiest way to explain an asymmetry in RM fluctuation amplitude that 
is correlated with jet sidedness is to postulate that the rotation arises in a 
large-scale gas component surrounding the source - most plausibly the group- 
scale hot component. The asymmetry is then simply due to the differin g path 
lengths through this gas to the approaching (brighter) and receding jets ( Laing! 



1988). In order to model the asymmetry we simulate the Faraday rotation from 



an isotropic, random magnetic field with a power spectrum corresponding to 
that derived from the RM structure-fu nction analysis emb edded in ionized gas 
with a smooth density distribution (cf. iMurgia et al. 20041 ). We first assume a 



spherically-symmetric density distribu tion with a core radius of 150 arcsec as fit 
to rosat data for the 3C31 group by iKomossa & Bohringer! ( 19991 ). An exam- 



ple RM distributio n for a radio source o f negligible thickness inclined by 52° 
to the line of sight ( Laing <fe Bridle! 120021 ) is shown in Figs^d) - (f). The rms 
central magnetic field strength for this model is (B 2 ) 1 ^ 2 = 0.21 nT (2.1 /zG). The 
predicted RM distribution shows an asymmetry, but there is no sharp change 
at the nucleus and the profile falls too rapidly with distance on the counter-jet 
side. The most plausible reason for the discrepancy is that the radio lobes have 
displaced the surrounding gas, which is therefore highly non-spherical. Cavi- 
ties in the X-ray-emitting gas asso ciated with radio lob es are indeed observed 
in similar sources such as 3C449 (|Croston et al.ll2003l 1. We assume that the 



contents of the cavity produce negligible Faraday rotation. If we model the 
cavity as initially conical, with a half-opening angle of ~55° within 100 arcsec 
of the nucleus and thereafter cylindrical, we can approximately reproduce the 
RM distribution (Figs B b - i). Such a cavity would not have been apparent in 



the existing X-ray data ( Komossa &: Bohringer 19991 ). but should be detectable 
with XMM-Newton. 
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